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A theory of angle-resolved photoemission (ARPES) in 
doped cuprates and other charge-transfer Mott insulators is 
developed taking into account the realistic (LDA+U) band 
structure, (bi)polaron formation due to the strong electron- 
phonon interaction, and a random field potential. In most 
of these materials the first band to be doped is the oxygen 
band inside the Mott-Hubbard gap. We derive the coher- 
ent part of the ARPES spectra with the oxygen hole spectral 
function calculated in the non-crossing (ladder) approxima- 
tion and with the exact spectral function of a one-dimensional 
hole in a random potential. Some unusual features of ARPES 
including the polarisation dependence and spectral shape 
in YBa2Cu307 and YBa2Cu40g are described without any 
Fermi-surface, large or small. The theory is compatible with 
the doping dependence of kinetic and thermodynamic prop- 
erties of cuprates as well as with the d-wave symmetry of the 
superconducting order parameter. 

PACS numbers:74.20.Mn,74.20.-z, 74.25. Jb 



1. Introduction 

The concept of polarons led to the discovery of the 
copper oxide superconductors. The expectation was that 
if 'an electron and a surrounding lattice distortion with 
a high effective mass can travel through the lattice as 
a whole, and a strong electron-phonon coupling exists, 
the perovskite insulator cauld be turned into a high tem- 
perature superconductor 'El. In this paper we develop the 
theory of ARPES in doped charge-transfer Mott insu- 
lators based on the bipolaron theoryBu which describes 
some unusual ARPES features of high-T c YBa2Cus07_5 
(Y123), YBa 2 Cu 4 8 (Y124) and a few other material^ 

ARPES data in cuprates remains highly controversial! 
One of the surprising features is a large Fermi surface 
claimed to exist in a wide range of doping fitting well the 
LDA band structures in the earlier studies. This should 
evolve with doping as (1 — x) in a clear contradiction 
with low frequency kinetics and thermodynamics (see, 
for exampleuil3, which show an evolution proportional 
to x io.-a.wide range of doping including the overdopcd 
regionoO (x is the number of holes introduced by dop- 
ing). Now it is established, however, that there is a nor- 
mal state (pseud<j4sap in ARPES and tunnelling, exist- 
ing well above T c E3ll3, so that some-segments of a 'large 
Fermi surface' are actually missingtfl. The temperature 
and doping dependence of the gap. still remain a sub- 
ject of controversy While kiatvtidj, thermodynamics, 
tunncllingE£l and some ARPESEJllj data suggest that the 
gap opens at any relevant temperature in a wide range 
of doping, other ARPES studiesEZrt£l claim that it ex- 
ists only in underdoped samples below a characteristic 



temperature T*. 

Perhaps the most intriguing feature of ARPES is 
the extremely narrow and intense peak lying below the 
Fermi energy, which most clearly seen near the Y 
and X points ia-Y124E], Y123EJ and more recently in 
La 2 _a;Sr a ;Cu04E3. Its angular dependence and spectral 
shape as well as the origin of the featureless (but dis=. 
persive) background remain unclear. Some authorsE2l 
refer to the peak as an extended van Hove singularity 
(evHs) arising from the plane (CUO2) strongly corre- 
lated band. They also implicate the resulting (quasi) ID 
density of states singularity as a possible origin for the 
high transition temperature. However, recent polarised 
ARPES studies of untwinned Y123 crystals of excep- 
tional qualitycj unambiguously refer the peak to as a nar- 
row resonance arising primarily from the quasi-lD CUO3 
chains in the buffer layers rather than from the planes. 
Also, careful analysis of the Eliashberg equations shows 
that the van Hove singularity can hardly be the origin 
of high T c , in sharp contrast with a naive weak-coupling 
estimate. Interestingly_a very similar narrow peak was 
observed by Park et aln in high-resolution ARPES near 
the gap edge of the cubic semiconductor FeSi with no 
Fermi surface at all. 

In this paper we take the view that cuprates and many 
other transition metal compounds are charge -transfer 
Mott-Hubbard insulators at any level of dopingcil. This 
means that the first band to be doped is the oxygen band 
lying within the Hubbard gap, Fig.l. Strong coupling 
with high frequency phonons, unambiguously established 
for many oxidesc3, leads to the high-energy spectral fea- 
tures of an oxygen hole in an energy window about twice 
the Franck-Condon (piglaronic) level shift, E p , and to the 
band-narrowing effecta. On the other hand, the low en- 
ergy spectral function is influenced by the low frequency 
thermal lattice, spin and random fluctuations. The lat- 
ter can be described as 'Gaussian white noise'. The p- 
hole polaron in oxides is almost one-dimensional due to 
a large difference in the ppa and ppn hopping integrals 
and effective 'one-dimensional' localisitaion by the ran- 
dom potential as described in Refrj. This allows us to 
explain the narrow peak in the ARPES spectra with the 
spectral density A(k,E) of a one-dimensional particle in 
a Gaussian white noise potentialcH. 

The electron-phonon interaction also binds holes 
into inteipsite oxygen bipolarons the size of a lattice 
constantoca. The bipolaron density remains relatively 
low (below 0.15 per cell) at any relevant level of dop- 
ing. The residual repulsive (Coulomb) polaron-bipolaron 
interactions is strongly suppressed by the lattice polar- 
isation owing to a large static dielectric constant. As 
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a result, the only role of hole bipolarons in ARPES is 
to pin the chemical potential inside the charge transfer 
gap, half the bipolaron binding energy above the oxygen 
band edge, Fig. la. This binding energy as well as the 
singlet-triplet bipolaron exchange energy are thought to 
be the origin of the japrmal state pseudogaps, as first 
proposed by one of usE3. In particular, the non-Korringa 
temperature dependence of NMRlj, and more recently 
the values and universal scaling -p.tJp_temperature of 
the uniform jaagnetic susceptibilityt^Jtrtil and electronic 
specific heated were quantitatively described with bipo- 
larons. In overdoped samples the bipolaron and polaron 
bands might overlap-because the bipolaron binding en- 
ergy becomes smallE^I, so the chemical potential might 
enter the oxygen band, Fig. lb, and then a Fermi- level 
crossing might be seen in ARPES. The featureless back- 
ground is explained as the phonon cloud of a small hole 
polaron, which spreads over a wide energy interval about 
2E p ~ leV. 

The same approach has been applied to the tunnelling 
spectraE]. As a result, the temperature independent 
gap and the anomalous gap/T c ratio, injection/emission 
asymmetry both in magnitude and shape, zero-bias con- 
ductance at zero temperature, the spectral shape inside 
and outside the gap region, temperature/doping depen- 
dence and dip-hump structure of the tunnelling conduc- 
tance were described. The approach is compatible with 
the recent |iLDA+U' band structures in cuprates and 
manganiteS, which suggest the single-particle density 
of states shown in Fig.l. It is clearly compatible with the 
doping evolution of thermodynamic and kinetic proper- 
ties because holes introduced by doping into the oxygen 
band are the only carriers in our theory. Moreover, the 
derived bipolaron energy dispersion with minima at the 
Brillouin zone boundaries presides d-w&ve symmetry of 
the Bose-Einstein cooeLensateES in agreement with phase- 
sensitive experimental. It also provides a parameter-free 
expression for T c in a few dozen cuprates, irrespective of 
the level of dopingL£lE3. 

Here we first derive the polaronic ARPES theory in 
Section 2. Then, in Section 3 we apply the simplest non- 
crossing approximation for the self-energy and compare 
it with the exact result in Section 4. The experimen- 
tal ARPES results, in particular the angular dispersion, 
spectral shape and polarisation dependence are compared 
with the theory in Section 5. We conclude that the nar- 
row peak in ARPES is an intrinsic polaron (oxygen) 
band found below the chemical potential by half of the 
bipolaron binding energy. We also arrive at the conclu- 
sion that the present experimental photocmission spec- 
tra can be understood with a small or even without any 
Fermi surface (depending on doping). 

2. Polaronic ARPES 

The interaction of the crystal with the electromagnetic 
field of frequency v is described by the Hamiltonian (in 
the dipole approximation) 



H int = {Zniy^sinivt) J^ie ■ dk,k')44' + H - c ; (!) 

k,k' 

where / is the intensity of the radiation with the po- 
larisation e (we take c=7i = feB = l), kis the mo- 
mentum of the final state (i.e. of the photoelectron 
registered by the detector), k' is the (quasi)momentum 
of the hole remaining in the sample after the emission, 
and and hy, are their creation operators, respec- 
tively. For simplicity we suppress the band index in h,, ■ 
Due to the translational symmetry of the Bloch states, 
k' >= u_k'(r)exp(— ik! ■ r), there is a momentum con- 
servation rule in the dipole matrix element, 

d k ,k' = d(k)<5 k+ k\G (2) 

with 

d(k)= le (A7« ) 1/2 V k I e- iGr u^ G (v)dr, (3) 

Jv 

where G is a reciprocal lattice vector, N the number of 
unit cells in the crystal of the volume Nvq. 

The Fermi Golden Rule gives the photocurrent to be 

I(k,E) =4vr 2 /|e-d(k)| 2 

x£e f2 +^-^|</ | hl_ G \i>\ 2 S(E + E f -E i ), (4) 

where E is the binding energy, Efj is the energy of 
the initial and final states of the crystal, and f2, /j,, Ni 
are the thermodynamic and chemical potentials and the 
number of holes, respectively. By definition the sum in 
Eq.(4) is n(E)A(k — G, —E) where the spectral function 
A(k -G,E) = (-l/7r)9G fl (k - G,E) is proportional to 
the imaginary part of the retarded Green function (GF), 
and n{E) — [exp(E/T) + the Fermi distribution. 

In the following we consider temperatures well below the 
experimental energy resolution, so that n{E) = 1 if E 
is negative and zero otherwise, and, for convenience, we 
put G = 0. 

The spectral function depends on essential interactions 
of a single hnlc with the rest of the system. As we ar- 
gued earlierc3 the most important interaction in oxides 
is the Frohlich long-range electron-phonon interaction of 
the oxygen hole with the c-axis polarised high-frequency 
phonons. The Frohlich interaction is integrated out with 
the Lang-Firsov canonical transfpcmationEj. With the 
Quantum-Monte-Carlo techniqueta one can prove that 
this transformation is practically exact for a long-range 
interaction in a wide region of the coupling strength, in- 
cluding the intermediate and weak-coupling regime and 
in a wide range of the adiabatic ratio, uj/D (lu is the 
phonon frequency and D is the bandwidth). By apply- 
ing the Lang-Firsov canonical transformation the hole 
Matsubara GF is expressed as a convolution of the coher- 
ent polaron GF and the multiphonon correlation function 
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where the multi-phonon correlation function, <r(m, r) = 
TJ2 n e _lf2 " T cr(m, Q n ) is found as 

a(m,T) = exp(±^\ 7 (q)\ 2 Um,T)\ . (6) 



Here 



/ q (m, r) = [cos(q • m)cosft,(w q |r|) - l]coth^ 
+ cos(q • m)sinh(ujq\T\) 

with m the lattice vector, uj n = irT(2n +1), n = 
0,±1,±2... and tt n = 2-rrTn. 

In the case of dispersionless phonons and a short-range 
(Holstein) interaction with a q -independent matrix ele- 
ment ((x> q = lj, |7(q)| 2 = 2g 2 ) one can readily calculate 
the Fourier component of er(m,T) to obtaind 
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The Green's function of a polaronic carrier comprises two 
different contributions. The first coherent k-dependent 
term arises from the polaron band tunneling. The spec- 
tral weight of the coherent part is strongly (exponen- 
tially) suppressed as Z = exp(—g 2 ) and the effective 
mass is strongly enhanced, £k = ZE^ — u (we include 
the polaronic level shift into the chemical potential, u) 
. Here E(k) is the bare (LDA+U) hole band disper- 
sion. The second k-independent contribution describes 
the excitations accompanied by the emission and absorp- 
tion of phonons. We believe that this term, Ii nco h{E) is 
responsible for the asymmetric background in the op- 
tical conductivity and in the photoemission spectra of 
cuprates and manganites. We notice that its spectral 
density spreads over a wide energy range of about twice 
the polaron level shift E p = g 2 ui. On the contrary the co- 
herent term shows an angular dependence in the energy 
range of the order of the polaron bandwidth 2w = ZD. 

As we have discussed above ARPES measures the 
imaginary part of the retarded GF multiplied by the 
Fermi-Dirac distribution function and by the square of 
the dipole matrix element. G R is obtained from Q(k, u> n ) 
with a substitution ioj n — » E + iO + . As a result we obtain 

I(k, E) ~ \d(k)\ 2 n(E)ZS(E + a) + hncohiE), (8) 

where Imcoh (E) is a structureless function of the binding 
energy, which spreads from about — u> down to —2E p . 
Only in the extreme limit of a very strong electron- 
phonon coupling, where the polaron bandwidth is well 
below the phonon frequency, can multiphonon structure 



in lincoh(E) be verified. Ji has actually been observed in 
the optical conduct ivityEj. 

The small Holstein polaron is very .heavy except if the 
phonon frequency is extremely higfc3. In fact, oxides 
are ionic semiconductors where the long-range electron- 
phonon interaction dominates. That leads, [to a much 
lighter mobile small Frohlich polaron (SFP)E3l3. Consid- 
ering the electron-phonon interaction in a multi-polaron 
system one has to take into account the dynamic prop- 
erties of the dielectric response functiorO. One can 
naively believec9 that the long-range Frohlich interaction 
becomes short range (Holstein) due to screening. This 
is of course untrue. In the long-wave limit (q — > 0) the 
response of polarons at the optical phonon frequency is 
dynamic, because u> » qv (v is the characteristic group 
velocity of polarons). Therefore, the singular (Frohlich) 
behaviour of 7(q) ~ 1/q is unaffected by the screening. 
Polarons are slow enough and cannot screen the high- 
frequency crystal field oscillations. Hence, the interaction 
with the high-frequency optical phonons in ionic polaron 
solids remains lops-range whatever the carrier density. 
It is easy to showed using Eq.(5) and Eq.(6) that for any 
finite-radius interaction with a (/-dependent matrix ele- 
ment the coherent part of ARPES spectra takes the same 
form as Eq.(8), but with a different spectral weight (Z) 
and effective mass (Z') renormalisation exponents. Also 
some k dependence of the incoherent background ap- 
pears if the matrix element of the electron-phonon inter- 
action depends on qcn. Hence, in general, 

I(k,E) ~ \d(k)\ 2 n(E)ZS(E + ^) +I incoh (k,E), (9) 

with the same Z — exp{—E p /uj) as in the case of the 
Holstein polaron but with the SFP bandwidth much less 
reduced, 6c = Z'E(k) — jj,, where Z' — exp(—jE p /u>). In 
general one finds 
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(10) 



and E p = (1/2) J2q l7(<j)| 2w ai, ^ ne numerical coefficient 
7 might be as small as 0.4E3 and even smaller in the 
cuprates with the nearest neighbour oxygen-oxygen dis- 
tance less than the lattice constant, 7 ~ 0.2t3, which 
is precisely confirmed by exact QMC calculation of the 
small polaron massEi. On the one hand this important 
result tells us that small polarons as well as intersite bipo- 
larons are perfectly mobile and can account for the high- 
T c values in cupratesEJ. On the other the coherent spec- 
tral weight remains strongly suppressed in polaronic con- 
ductors, Eq.(9), because Z might be less than Z' by one 
or even a few orders of magnitude. These unusual SFP 
spectral features provide an explanation for the apparent 
discrepancy between a very small Drude weight and a rel- 
atively moderate mass enhancement, m* ~ 3m (r ^ri-J.0m e 
(depending on doping) of carriers in manganitest3'E-3 and 
cuprates— Hiey also explain why the evHs observed in 
ARPESl3l3 can be hardly seen in angle averaged pho- 
toemission. Indeed, the integrated spectral weight of the 
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incoherent background is proportional to (1 — Z), Eq.(7). 
It turns out to be much larger than the coherent contri- 
bution, proportional to Z « 1. Finally, the k depen- 
dent incoherent background might obscure the experi- 
mental determination of the Fermi-lepiel crossing, leaving 
the 'leading edge' determination of itlia rather unreliable. 

In the following we concentrate on the angular, spec- 
tral and polarisation dependence of the first coherent 
term, Eq.(9). The present experimental resolution.] al- 
lows probing of the intrinsic damping of the coherent 
polaron tunnelling. This damping appears due to the 
random field and low-frequency lattice and spin fluctu- 
ations described by the polaron self-energy S p (k, E), so 
that the coherent part of the spectral function is given 
by 



A p (k,E) 



-(1/tt) 



3E p (k, E) 



[E + m p (k,E)-^)Y 



[3£ p (k,£)F 
(11) 



Hence, the theory of the narrow ARPES peak reduces 
to the determination of the self-energy of the coherent 
small hole polaron scattered by impurities, low- frequency 
deformation and spin fluctuations. 

3. Self-energy of ID hole in the 
non-crossing approximation 
Due to energy conservation small polarons exist in 
the Bloch states at temperatures well below the optical 
phonon frequency T << to/ 2 no, matter how strong their 
interaction with phonons isEllc3cZln. This textbook re- 
sult, known for a long time, has been questioned by some 
authorso. It has been showrO that the confusion is due 
to a profound misunderstanding of the strong-coupling 
expansion by those authors. The finite polaron self- 
energy appears only due to the (quasi)elastic scattering. 
First we apply the simplest non-crossing (ladder) approx- 
imation to derive the analytical results, Fig. 2. Within 
this approximation the self-energy is k -independent for 
a short-range scattering potential like the deformation or 
a screened impurity potential, so that 



£ p (£)~]TG«(k,£), 



(12) 



with G*(k, E) = [E-&- ^p(E)}- 1 

The oxygen polaron-spectrum is parametrised in the 
tight-binding model ast£l 



f 1 - 



2[tcos{k x>y a) — t'cos(k y , x a) + t c cos(k z d)] — fi, 



(13) 



If the oxygen hopping integrals in Eq.(13), reduced by 
the narrowing effect, are positive, the minima of the po- 
laron bands are found at the Brillouin zone boundary 
in X (%, 0) and Y (0, tt). The wave vectors correspond- 
ing to the energy minima belong to the stars with two 
prongs. Their group has only ID representations. This 
means that the spectrum is degenerate with respect to 



the number of prongs of the star. The spectrum Eq.(13) 
belongs to the star .-with two prongs, and, hence it is a 
two-fold degenerateLB. The doublet is degenerate if the 
hole resides on the apical oxyger£3. In general, the de- 
generacy can be removed due to the chains in the buffer 
layers of Y123 and Y124, so that the y-polaronic band 
corresponding to the tunneling along the chains might be 
the lowest one. 

As mentioned above the oxygen hole is (quasi) one- 
dimensional due to a large difference between the oxy- 
gen hopping integrals for the orbitals elongated paral- 
lel to and perpendicular to the oxygen-oxygen hopping 
t',t c << t. This allows us to apply a one-dimensional 
approximation, reducing Eq.(13) to two ID parabolic 
bands near the X and Y points, = k 2 /2m* — /x with 
to* = l/2ia 2 and k taking relative to (71", 0) and (0,7r), 
respectively. Then, the equation for the self-energy in the 
non-crossing approximation, Eq.(12) takes the following 
form 



Sp(e) 



-3/2 



P P (e) 



-1/2 



(14) 



for each doublet component. Here we introduce a dimen- 
sionless energy (and self-energy), e = (E + u)/eo using 
€q = (D 2 m*) 1 ^ 3 as the energy unit. The constant D 
is the second moment of the Gaussian white noise po- 
tential, comprising thermal and random fluctuations as 
D = 2(VqT/M + n lm v 2 ), where V is the amplitude of 
the deformation potential, M is the elastic modulus, n,- m 
is the impurity density, and v is the coefficient of the 6- 
function impurity potential. The solution is 



S P (c) 



1 



i&l 2 



l-iZ 1 ' 2 



1 

16 



1 

16 4 

£l 

" 27 



e 

27 



(-L 

V 256 



1 

256 



e 

216 



216 

1/2-1 V3 



1/3 



(15) 



While the energy resolution in the present ARPES 
studies is almost perfectta, the momentum resolution re- 
mains rather poor, 6 > O.ln/a. Hence we have to in- 
tegrate the spectral function, Eq.(ll), with a Gaussian 
momentum resolution to obtain the experimental pho- 
tocurrent, 



I(k,E) - Z dk'A p (k',-E)exp[- 



(fc - k ') 



A2 



(16) 



The integral is expressed in terms of S p (e), Eq.(15) and 
the tabulated Error function w(z) as 



I(k,E) 



2Z, 



C£ p (e)[w( Zl )+w(z 2 )}), (17) 



where Z1.2 = [ik — i/2£ p (e)]/<5, w(z) = e~ z erfc(-iz) 
and e = (— E + u)/e.Q. This photocurrent is plotted as 
dashed lines in Fig. 3 for two momenta, k = 0.047r/a (al- 
most Y or X points of the Brillouin zone) and k = O.Sn/a. 
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The chemical potential is placed in the charge transfer 
gap below the bottom of the hole band, fi = —20meV, 
the momentum resolution is taken as 5 — 0.287r/a and 
the damping eo — 19rneV. 

The imaginary part of the self-energy, Eq.(15) disap- 
pears below e = — 3/2 5 / 3 ~ —0.9449. Hence this approxi- 
mation gives a well-defined gap rather than a pseudogap. 
Actually, the non-crossing approximation fails to describe 
the localised states inside the gap (i.e. the Lifshitz tail of 
the density of states). One has to go beyond the simple 
ladder, Fig. 2, .to describe the single-electron tunnelling 
inside the gapt2l and the ARPES spectra at very small 
binding energy. 

4. Exact spectral function of ID hole 

The exact spectral function for a one dimensional par- 
ticle in a random Gaussian white noise potential was de- 
rived by rlaJperinEZI and the density of states by Frisch 
and Lloydcil. Halperin derived two pairs of differential 
equations from whose solutions the spectral function of 
a 'Schrodinger' particle (i.e. in the effective mass ap- 
proximation) and of a 'discrete' particle (tight-binding 
approximation) may be calculated. Thep-QMC polaronic 
bandwidth is about 100 meV or largerEa, which allows 
us to apply the 'Schrodinger' particle spectral function, 
given byEZl 



Ap(k, e) = 4 / p (-z)$lpi(z)dz, 



where po.i(z) obey the two differential equations 

Po = 0, 



d 2 d , 2 „ ; 



and 



d 2 d . , „ . 
-t— 2 + — (z 2 + 2e) - z - ik 



dz 2 ' dz ' 
with boundary conditions 



pi +p = 0, 



lim z 2 n p n (z) = lim z 2 n p n {z) 



(18) 



(19) 



(20) 



(21) 



where k is measured in units of ko = (D^^m*) 2 ! 3 . The 
first equation may be integrated to give 



Po(z) 



exp{-z 3 /3 - 2ze) /'^ exp(u 3 /3 + 2ue)dt 
tt 1 / 2 /°° u" 1 1 12 exp(-u 3 / '12 - 2ue)du 



(22) 



The equation for p\(z) has no known analytic solution, 
and hence must be solved numerically. There is however 
an asymptotic expression for A p (k, e) in the tail where 
|e|»l: 



A p (k,e) ~ 27r(-2e)5 exp[--(-2e)^]cosh 2 



irk 



(-8e)3 



(23) 



In practice, for computational efficiency we use the exact 
spectral density for — 1.4 < e < 1, and outside this range 
we use the asymptotic result for e < —1.4, Eq.(23) and 
the non-crossing approximation for e > 1, where they 
are almost indistinguishable from the exact result on the 
scale of the diagrams plotted here. 

The result for A p (k, —E) integrated with the Gaussian 
momentum resolution is shown in Fig. 3 for two values of 
the momentum (solid lines). Quite differently from the 
non-crossing approximation the exact spectral function 
(averaged with the momentum resolution function) has 
the Lifshitz tail due to the states localised by disorder 
within the normal state gap. However, besides this tail 
the non-crossing approximation gives very good agree- 
ment, and for binding energy greater than about 30meV 
it is practically exact. 

The cumulative DOS 



N p (e) = (2tt)- 



de' 



dkA p (k,e') (24) 



is expressed analyticallyEd in terms of the tabulated Airy 
functions Ai(x) and Bi(x) as 



N p (e)) = tt- 2 [Ai 2 {-2e) + Bi 2 {-2e)} 



(25) 



The DOS dN(e)/de fits well the, voltage-current tun- 
nelling characteristics of cuprates£3. 

5. Theory of ARPES in Y124 and Y123 

With the polaronic doublet, Eq.(13) placed above the 
chemical potential we can auantitativehi, describe high- 
resolution ARPES in Y123B and Y124E3. First we ex- 
plain the experimentally observed polarisation depen- 
dence of the intensity near Y and XEil. The Bloch pe- 
riodic function Mk(r) can be expressed in terms of the 
Wannier orbitals w(r) as 



u k (r) = N- 1/2 J2 e lk ' (m ~ r Mr - m). 



(26) 



Then the dipole matrix element is given by the deriva- 
tive of the Fourier component of the atomic (Wannier) 

orbital, /k = v^ 1 ^ 2 J dvw(r)exp(ik ■ r) as 



d(k) = i(e • V k )/ k . 



(27) 



To estimate /k we approximate the x, y oxygen orbitals 
contributing to the x and y polaronic bands, respec- 
tively, with w x (r) — (l/8alTr) 1 / 2 (x/2a )exp(-r/2a ) 
and w y (r) — [1 / Sa^ir) 1 / 2 (y /2ao)exp(— r/2ao). As a re- 
sult we obtain for the x orbital, 

|^ = (8a^/« ) 1/2 a 

X x ([(fca ) 2 + l/4]- 3 -6(A :;ca o) 2 [(A : ao) 2 + l/4]- 4 ), (28) 



and 



" Jk =-^alix/v ) 1 ' 2 alk x k y>z [{ka f + l/A}- i , (29) 



Ok 
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Here k is the photoelectron momentum and ao is the size 
of the Wannier function. For the case of y-orbital one 
should interchange x and y. Near the X and Y points 
of the Brillouin zone, \k y _ x \ << k, respectively. Then it 
follows from Eq.(28) andEq.(29) that the ARPES peak 
should be seen at X and almost disappear at Y if the 
photons are polarised along the x-direction, i.e. e || a. 
If the polarisation is along the y-direction (e || b) the 
peak appears at Y and almost disappears at X. Precisely 
this behaviour is observed in ARPES spectra obtained 
using polarised photonscll, Fig. 4. We also notice a very 
strong dependence of the dipole matrix element, Eq.(28) 



on the photon energy, d 



at large v. Hence, it is 



not surprising if the ARPES pjak disappears at large v 
as has been recently observedEj. 

The exact ID polaron spectral function, Eq.(18), inte- 
grated with the experimental momentum resolution, pro- 
vides a quantitative fit to the ARPES spectra in Y124 
along the Y — V direction, as shown in Fig. 5. The 
angular dispersion is described with the polaron mass 
to* = 9.9to e in agreement.. with the Monte-Carlo calcu- 
lations of the SFP massEa. The spectral shape is re- 
produced well with eo = 19meV, Fig. 6, in close agree- 
ment with tbe value of this parameter found in tunnelling 
experimental. That yields an estimate of the polaron 
scattering rate, which appears to be smaller than the po- 
laron bandwidth-/ about lOOmeV or larger), in agreement 
with the notionEj that many high-T c cuprates are in the 
clean limit. There is also quantitative agreement in the 
perpendicular direction Y — S, Fig. 7, in a restricted re- 
gion of small k x , where almost no dispersion is observed 
around Y. Slight dispersion in the Y — S direction to- 
wards the chemical potential might be due to a negative 
t' in Eq.(13). 

However, there is a significant loss of the energy- 
integrated intensity along both directions, Fig. 8, which 
the theoretical spectral function alone cannot account 
for. The energy- integrated ARPES spectra obey the sum 
rule, 



dEI(k,E) ~ |d(k)| 2 n k , 



(30) 



where n^ = (h^h^). If the dipole matrix element is al- 
most k-independent and the chemical potential is pinned 
well inside the charge-transfer gap, so that rik = 1 this 
integral would be k independent as well. This is not the 
case for Y124, no matter what the scanning direction 
is, Fig. 8. Therefore, we have to conclude that either the 
dipole matrix element is k dependent or (and) the oxygen 
band is strongly correlated (in the Mott-Hubbard sense) . 
As we have mentioned above, the incoherent background 
of SFP is angle dependent as well, which might contribute 
to the intensity loss. 

The rapid loss of the integrated intensity in the Y — Ji 
direction was interpreted by Randeria and Campuzancc3 
as a Fermi-surface crossing. While a Fermi-surface cross- 
ing is not incompatible with our scenario (see Fig.l in- 



set), we do not believe that it has really been observed 
in Y124. First of all these authors suppressed a few ex- 
perimental curves in the Y — T direction, Fig. 5, which 
prevented them from observing the intensity loss in this 
'dielectric' direction, where there is obviously no Fermi- 
surface crossing. This loss of intensity along Y — T tells 
us that the intensity loss might be due to the matrix el- 
ement rather than to the Fermi-surface crossing in both 
directions. This is confirmed by our observation of_a 
similar rapid loss of the intensity in a dielectric, FeSiEij, 
Fig. 9 with no Fermi-surface at all. The peaks in the 
Y — S direction are all 15 meV or more below the Fermi 
level - at a temperature of 1 meV, if the loss of spectral 
weight were due to a Fermi-surface crossing one would 
expect the peaks to approach much closer to the Fermi 
level. Also the experimental spectral shape of the in- 
tensity at k = kp is incompatible with any theoretical 
scenario, including different marginal Fermi-liquid mod- 
els, Fig. 10. The spectral function on the Fermi surface 
should be close to a simple Lorentzian, 

\E\P 



A p (k F ,E) 



E 2 + constant x E 2/3 ' 



(31) 



because the imaginary part of the self-energy behaves as 
\E\" with < (3 < 2. On the contrary, the experimental 
intensity shows a pronounced minimum at the alleged 
Fermi-surface, Fig. 10. 

If there is indeed no Fermi-surface crossing in many 
cuprates, as we argue, why then does the 'maximum 
locus' determination point to a large Fermi surface in 
cuprates, which is drastically incompatible with their ki- 
netic and thermodynamic properties? We propose that it 
appears due to the fact that oxygen semiconducting band 
has its minima at large k inside or even on the boundary 
of the Brillouin zone. That is why ARPES show intense 
peaks near large k imitating a large Fermi-surface. 
6. Summary and conclusions 

In summary, we have proposed a theory of ARPES 
in cuprates based on the LDA+U band structure and 
the bipolaron theory compatible with the normal state 
kinetic and thermodynamic properties of these materi- 
als. The theory explains the narrow flat bands observed 
in Y123 and Y124, including their polarisation, spectral 
and angular dependence, as well as a featureless (but dis- 
persive) background. The ARPES peak originates from 
the hole excitations of the polaronic oxygen band of the 
buffer layers, in agreement with the experimental results 
and electronic structure of Schabel et aZcll. Differently 
from these authors we suggest that this band is intrin- 
sic for cuprates and takes part in the bipolaron formation 
and superconductive, which is nicely confirmed by a few 
independent studies^!! The normal state gap is half of 
the bipolaron binding energy. The angular dependence 
of the peak and of the gap is due to the polaron band 
dispersion, which agrees well with the QMC results for 
the small Frohlich polaron. 

The spectral shape of the peak is affected by the soft 
lattice, spin and random fluctuations. The characteristic 



G 



scattering rate agrees well with that found in the tun- 
nelling experimental. This scattering rate is tempera- 
ture dependent not only due to the thermal lattice fluctu- 
ations, but also because of the anomalous screening below 
T c in the charged Bose-liquida. The Bose-Einstein con- 
densate screens effectively the long-range Coulomb po- 
tential of impurities. As a result one can expect a drastic 
change of the damping eo when T c is passed. That might 
help to understand the near disappearance of the narrow 
peak above T c in some Bi-cuprates. On the other hand in 
the stehiometric Y124 with (theoretically) no impurities, 
one can expect about the same results from ARPES be- 
low and above T c , which seems to be the case£3. Within 
the bipolaron theory there is only one single — particle 
gap, which is half of the bipolaron binding energy both 
below and above T c . 

We believe that many cuprates are doped insulators 
with no Fermi surface at all due to the bipolaron forma- 
tion. The Fermi-surface crossing, if it were firmly estab- 
lished in the overdoped samples, would correspond to a 
small Fermi surface of the oxygen band pockets located 
at finite k like in many ordinary semiconductors, for ex- 
ample, in Ge and Si. 

The authors greatly appreciate enlightening discus- 
sions with D.S. Dessau, N.E. Hussey, V.V. Kabanov, 
P.E. Kornilovitch, G.J. Kaye, A.I. Lichtenstein, G.A. 
Sawatzky, J.R. Schrieffer, Z.-X. Shen, J. Zaanen, G. 
Zhao, and R. Zeyher. C.J.D. has been supported in this 
work by a grant from the EPSRC of the UK. 
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Figure Captures 

Fig.l. Schematic LDU+U density of states. The chem- 
ical potential is pinned inside the charge transfer gap (a) 
due to the bipolaron formation in underdoped cuprates. 
It might enter the oxygen band in overdoped cuprates 
(b) if the polaron band crosses the bipolaron one (inset). 

Fig. 2. The non-crossing diagram for the self-energy. 
The dashed line corresponds to the random potential and 
(or) to the thermal lattice and spin fluctuations. 

Fig. 3. The polaron spectral function, integrated with 
the momentum resolution function for two angles, k = 
O.OAir/a (upper curves), and k = 0.307r/a with the damp- 
ing eo = 19meV, the momentum resolution 6 = 0.287r/a 
and the polaron mass m* = 9.9m e . The bipolaron bind- 
ing energy 2\n\ = AOmeV. The dashed curves are the 
spectral density integrated with the momentum resoltion 
in the non-crossing approximation. 

Fig-4. Polarisation dependence of the ARPES peak in 
Y1230 near X and Y points. 

Fig. 5. Theoretical ARPES spectra (b) compared with 
experiment (a) in Y124c2lfor Y — T direction. Parameters 
are those of Fig. 3. 

Fig. 6. Theoretical fit (dashed lines) to two experimen- 
tal ARPES curves corresponding to k — 0.047r/a (upper 
curves), and k = 0.307r/a. 

Fig. 7. Theoretical ARPES spectra (b) compared with 
experiment (a) in Y124E9 for Y — S direction. 

Fig. 8. The energy-integrated ARPES intensity in Y124 
in the Y — T (a) and Y — S (b) directions. Momenta are 
measured relative to the Y-point of the Brillouin zone. 

Fig,9. (a) The energy-integrated ARPES intensity in 
FeSicJ as a function of the analyser angle. The spectra 
are shown in (b). 

Fig. 10. The experimental ARPES signal (solid line) 
on the alleged Fermi-surface does not correspond to a 
Fermi- liquid spectral function (dashed line). We assume 
particle-hole symmetry to obtain the spectral function 
for negative binding energy. 
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